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Abstract 

The big bounce (BB) transition within a flat Friedmann-Robertson- Walker model is analyzed in 
the setting of loop geometry underlying the loop cosmology. We solve the constraint of the theory at 
the classical level to identify physical phase space and find the Lie algebra of the Dirac observables. 
We express energy density of matter and geometrical functions in terms of the observables. It is 
the modification of classical theory by the loop geometry that is responsible for BB. The classical 
energy scale specific to BB depends on a parameter that should be fixed either by cosmological 
data or determined theoretically at quantum level, otherwise the energy scale stays unknown. 
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I. INTRODUCTION 



It is commonly believed that the cosmological singularity problem [j], 0, [3| may be resolved 
in a theory which unifies gravity and quantum physics. Recent analyses done within the loop 
quantum cosmology (LQC) concerning homogeneous isotropic universes of the Friedmann- 
Robert son- Walker (FRW) type, strongly suggest that the evolution of these universes does 
not suffer from the classical singularity Strong quantum effects at the Planck scale cause 
that classical big bang is replaced by quantum big bounce (BB) 0, H, @, 0, 0]. 

The resolution of the cosmic singularity problem offered by LQC requires the existence 
of a fundamental length, which effectively implies the discreteness of quantum geometry. 
However, the size of this length has not been determined satisfactory yet, i.e. derived within 
LQC. Presently, it is an ad-hoc assumption of standard LQC 0, Ho| . 

Our paper is an extended version of the classical part of [ljj . Its goal is the demonstration 
that the resolution of the initial cosmological singularity is due to the modification of the 
classical theory by the loop geometry. Quantum effects seem to be of a secondary importance, 
but should be examined since the energy scale specific to BB has not been identified yet 

llj]. The big bounce may occur deeply inside the Planck scale where it is commonly 
expected that quantum effects cannot be ignored. 



The difference between standard LQC [LJ, [15| and our nonstandard LQC is the following: 



(i) we determine an algebra of observables on the kinematical phase space; (ii) we solve 
the Hamiltonian constraints at the classical level to identify physical phase space (i.e. the 
space of Dirac's observables); (iii) we express functions on the physical phase space, like 
the matter density and geometrical operators (length, area, volume), in terms of Dirac's 
observables and an evolution parameter; and (iv) by quantization we mean: finding a self- 
adjoint representation of the algebra of the Dirac observables and solution to the eigenvalue 
problem for operators corresponding to functions specified in (iii). Roughly speaking, in 
standard LQC one does not identify an algebra of physical observables and one imposes the 
Hamiltonian constraints only at the quantum level. 

The present paper concerns the classical level so we are mainly concern with items (i)- 



(iii). Our next paper [16j| will be an extended version of the quantum part of [111 ], i.e. it will 
be devoted to the realization of item (iv). 

For simplicity of exposition we restrict ourselves to the flat FRW model with massless 
scalar field. This model of the universe includes the initial cosmological singularity and has 
been intensively studied recently within LQC. 

In order to have our paper self-contained, we recall in Sec. II the form of classical 
Hamiltonian in terms of holonomy-flux variables. In Sec. Ill we solve the constraint and 
analyze the relative dynamics in the physical phase space. Section IV is devoted to the 
algebra of observables. We consider the energy density of the scalar field and geometrical 
operators in Sec. V. Higher order holonomy corrections are shortly discussed in the appendix. 
We conclude in the last section. 
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II. HAMILTON! AN 



The gravitational part of the classical Hamiltonian, H gi in gen eral relativity is a linear 
combination of the first-class constraints, and reads [13, EH, 23, 24, 2^ 

H g := [ d 3 x(N l d + N a C a + NC), (1) 

where £ is the spacelike part of spacetime RxE, (N l , N a , N) denote Lagrange multipliers, 
(Ci,C a ,C) are the Gauss, diffeomorphism and scalar constraint functions. In our notation 
(a, b = 1, 2, 3) are spatial and (i, j, k = 1, 2, 3) internal SU{2) indices. The constraints must 
satisfy a specific algebra. 

Having fixed local gauge and diffeomorphism freedom we can rewrite the gravitational 
part of the classical Hamiltonian (for the flat FRW model with massless scalar field) in the 
form (see, e.g. [7|) 

H g = -r 2 [ d 3 x Ne-'e^E^F^ , (2) 



JV 

where 7 is the Barbero-Immirzi parameter, V C X is an elementary cell, X is spacelike 
hyper-surface, N denotes the lapse function, e^ k is the alternating tensor, Ef is a densitized 
vector field, e := \J | det E | , and where F* b is the curvature of an SU(2) connection A l a . 

The resolution of the singularity, obtained within LQC, is based on rewriting the curvature 
F^ b in terms of holonomies around loops. The curvature F^ b may be determined [3] by making 
use of the formula (see the appendix) 

fjM _ 1 

where 

h^ = ^W)- 1 (h^)- 1 (4) 
is the holonomy of the gravitational connection around the square loop Dy, considered over 

-1/0 

a face of the elementary cell, each of whose sides has length \iV (where /i > 0) with 
respect to the flat fiducial metric °q a b '■= % °uj l a ijj^ a \ fiducial triad °e\. and cotriad °ui k satisfy 
°ujI "e'j = Sj] the spatial part of the FRW metric is q a b = a 2 (t) °q ab ; Ar denotes the area 



of the square; V a = j v ^f°qd?x is the fiducial volume of V. Figure 1 shows geometrical setup 
for determination of . 

LJjj 

1 /3 

The holonomy along straight edge °e1d a of length fiV reads 

h^\c) = Vexp ( / r {k) A ( ^dx a ) = exp(r fc /ic) = cos(/xc/2) 1 + 2 sin(//c/2) r k , (5) 
Jo 



where r k = —ia^jl (a k are the Pauli spin matrices) and V denotes the path ordering symbol. 
Equation (5) presents the holonomy calculated in the fundamental, j = 1/2, representation 
of SU(2). The connection A\ and the density weighted triad E% which occurs in (jSJ) is 
determined by the conjugate variables c and p as follows: A\ = °u*cV ~ 1/3 and E^ = 

° e t y/QoPVo 2//3 ; where c = 'jdV^ 3 and \p\ = a 2 V^ 3 . Equation (jSj) presents the holonomy 
calculated in the fundamental, j = 1/2, representation of SU(2). 
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FIG. 1: Holonomy of connection around the square loop. Suitable fiducial cotriads °ui^ are chosen 
to diagonalise the connection A*. 



Making use of (J2J), (J3J) and the so-called Thiemann identity [2 

e ijk e^E^E bk = ^ £ ^^Tr(^{(^)- 1 ,n r,) (6) 
27rG'7/iy / k v ' 

leads to H g in the form 



^ = lim ffW (7) 
where 

^ = -^0^Y, N ^ ( 8 ) 

and where V = = a 3 Vo is the volume of the elementary cell V. 

The classical total Hamiltonian for FRW universe with a massless scalar field, 0, reads 

H = H g + H (j> ^0, (9) 

where H g is defined by (J7J). The Hamiltonian of the scalar field is known to be: = 
N p^\p\~2 /2, where and p^ are the elementary variables satisfying {(j),p^\ = 1. The 
relation H m defines the physical phase space of considered gravitational system with 
constraints. 

Making use of (jSJ) we calculate © and get the modified total Hamiltonian ifg^ corre- 
sponding to (Q in the form 



where 



^ : =p7^ v: =IpI 3/2 (ii) 



are the canonical variables proposed in [3]. The variable (3 = 'yd/ a so it corresponds to 
the Hubble parameter a/a, whereas v 1 ^ 3 = aV^ 3 is proportional to the scale factor a. The 
relationship between the coordinate length \i (which depends on p) and the physical length 
A (which is a constant) reads 

(12) 



A = [i Ip] 1 ^ 2 = /id Vq 



1/3 



The complete Poisson bracket for the canonical variables (f3,v, 4>iV^) is defined to be 

d- d- d- d- 



{•,-}:=47rG 7 



8(3 dv 



dv d/3 



+ 



dptf, dptf, 



The dynamics of a canonical variable £ is defined by 
where £ 



(13) 



(14) 



d^/dr, and where r is an evolution parameter. The dynamics in the physical 
phase space, J-p^ ys , is defined by solutions to f fT4"l) satisfying the condition pa 0. The 

solutions of (THI) ignoring the constraint pa are in the kinematical phase space, . Fffl . 

In what follows we apply the Dirac method of dealing with Hamiltonian constraints [17| : 
the Poisson bracket is worked out before one makes use of the constraint equations. 



III. DYNAMICS 

In the case a Hamiltonian is a constraint which may be rewritten in the form of a product 
of a simpler constraint and a function on J-j^ which has no zeros, the original dynamics 
may be reduced (to some extent) to the dynamics with the simpler constraint. 

Equation (TTU]) can be rewritten as 

H w = NH ( X) H W pa 0, (15) 

where 

where k 2 = 47rG/3. 

It is clear that = only in the case when p^ = = sin(A/3). Such case, due to (l2"Tj) - ([2"5l) . 
implies no dynamics. 

A. Relative dynamics 

An equation of motion for a function / defined on physical phase space, due to (fl4"|) . reads 
/ = {f,NH^H^} = {f,NH { X) }H^ + NH { X) {f,H^} = NH { Q X) {/, (17) 
since = 0. By analogy, for other function g we have 

g = {g, NH { X) H^} = NH { Q X) {g, H^}, for if W Pa 0. (18) 
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Therefore we have the relation 



which we rewrite in the form 

df dg 



(20) 

{f,HW} {g,HW} 

Equation (|20|) shows that in the case of the relative dynamics we have: (i) one phase space 
variable may be used as an 'evolution parameter' of all other variables, (ii) dynamics is 
gauge independent in the sense that there is no dependance on the specific choice of the 
lapse function N, and (iii) suitable choice of N may lead to a simpler form of Hamiltonian. 



B. Solution of the relative dynamics 

Since the relative dynamics is gauge independent, it is reasonable to chose the gauge 
which simplifies the calculations. Our choice is N := 1/H^ X \ In this gauge the equations of 
motion read 

P* = 0, (21) 

p = -4.G 7 (22) 

= K7 sgn(p^), (23) 
v = AnG-fv cos(A (3) sgn(sin(A /?)), (24) 
H {x) = 0. (25) 

Combining (123]) with (124"]) gives 

v 

— = 3kv cos (A/?) sgn(sin (A/?)) sgn(p0). (26) 

<\> 

Rewriting (126]) (and using v/<j>= dv/dcj)) gives 

sgn(sin(A/5)) dv 

__ = 3k sgn(p ) d$ 27 

cos(Ap) v 

Making use of the identity sin 2 (A/?) + cos 2 (A/3) = 1 and ( 125]) gives 



cos 



m = \/l-(^) 2 (28) 



Combining (T27]) with ([28]), for (3 e]0, tt/2A[, leads to 

dv 



3k sgn(p^) (29) 



a/v 2 - (^Ap^) 2 

Since is just a constant (due to (121"]) ) we can easily integrate (129]) and get 



In 



v + \Jv 2 - (n-fXp^y 



3k sgn(p )(0 - fa). (30) 
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Rewriting (130]) leads to 

2v = exp (3k sgnQ^) - O )) + (^M^) 2 • exp ( - 3k sgn(p ) (0 - O ))- (31) 
The solution for the variable (3 may be easily determined from (125]) rewritten as 



|sin(A/3)| , . 

^Tl^l = ^ J (32) 



The final expression reads 



sin(A/3) = — , 2^AN (33) 

exp (3Ksgn(p^) (0 - O )J + (ftr/M^exp ( - 3Ksgn(p^) (4> - O )) 

where the domain of the variable f3 has been extended to the interval ]0, 7r/A[. 

Equations fl3Tl) and fl33l) present the dependence of the canonical variables v and j3 on 
the evolution parameter <f), which is a monotonic function due to (123]) . 



IV. ALGEBRA OF OBSERVABLES 

A function, O, defined on phase space is a Dirac observable if 

{O,H (x) }^0. (34) 

Since we have 

{O, # (A) } = {O, NH^HW} = NH { x) {O, H w } + {£>, NH^}H^ X \ (35) 

it is clear that on the constraint surface, = 0, the Dirac observable satisfies (indepen- 
dently on the choice of N) a much simpler equation 

{O,H (x) }^0. (36) 
Thus, we put N := and solve (134]) in the whole phase space, i.e. we solve the equation 

— - v cos(A/3) — — — = 0. (37) 

A function O = 0(O u ...O h ) satisfies §H) if 

{Ox, ffW} = = {0 2 , F«} = . . . = {O k , H™}, (38) 
where k + 1 is the dimension of the kinematical phase space. It is so because one has 

{O, m = H {0 lt HW) + . . . + ^ {O fc , (39) 

In what follows we consider only elementary observables. The set of such observables, £ , 
is defined by the requirements: (i) each element of £ is a solution to (137]) . (ii) elements of £ 
are functionally independent on the constraint surface, = 0, (iii) elements of £ satisfy 
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a Lie algebra, and (iv) two sets of observables satisfying two algebras are considered to be 
the same if these algebras are isomorphic. 

In our case k = 3 and solutions to (1371) are found to be 



0i:=2V, 2 :=0-^arth(cos(A/3)), 3 :=sv^^-, (40) 

where s := sgn(p^). One may verify that the observables satisfy the Lie algebra 

{0 2i O x } = l, {O 1 ,O 3 } = 0, {0 2 ,0 3 } = 1 k. (41) 
Because of the constraint = (see (1321) ). we have 

3 = 7«0i- (42) 

Thus, we have only two elementary Dirac observables which may be used to parameterize 
the physical phase space T^ ys . To identify the Poisson bracket in J\^ yB consistent with the 

Poisson bracket (!T3"j) defined in J 7 ^, we find a symplectic twoform corresponding to (TT31 . 
It reads 

lu = — — <i/3 A dv + d(f) A rfp^. (43) 

47TG7 

The twoform u; is degenerate on due to the constraint = 0. Making use of the 

explicit form of this constraint (1321) and the functional form of 0\ and 0%, leads to the 
symplectic form Q on J r p'^ ys - Direct calculations give (see App. B) 

n-.= LU imx)=0 = dO 2 AdO u (44) 

where LUijj(\) =0 denotes the reduction of lu to the constraint surface. The Poisson bracket 
corresponding to ( 1441) reads 

d- d- d- d- 

so the algebra satisfied by 0\ and 2 has a simple form given by 

{0 2 ,0 1 } = l. (46) 

Our kinematical phase space, J 7 ^, is four dimensional. In relative dynamics one variable 
is used to parametrize three others. Since the constraint relates the variables, we have 
only two independent variables. This is the reason we have only two elementary physical 
observables parametrizing J\^ ya - 



V. FUNCTIONS ON PHASE SPACE 

In this section we discuss the functions on the constraint surface that may describe singu- 
larity aspects of our cosmological model. Considered functions are not observables, but they 
can be expressed in terms of observables and an evolution parameter 0. They do become 
observables for each fixed value of 0, since in such case they are only functions of observables. 
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A. Energy density 



An expression for the energy density p of the scalar field <fi reads 

2 

In terms of elementary observables we have 

p 4> = 1 , v = K7A|Oi| cosh(3/c(0- 2 )). (48) 

For fixed p^ the density p takes its maximum value at the minimum value of v. Rewriting 
( ]3Tj) in the form 

v 

— = cosh (3/ts(0 — 0o) — In A), where A := K7A (49) 

we can see that cosh(-) takes minimum value equal to one at 3ks (<j) — 4>o) = In A. Thus, 
the maximum value of the density, p max , corresponds to v — A and reads 

Pmax = 9 2 "\2" (^0) 

We can determine p max if we know A. However, A is a free parameter of the formalism. 
Thus, finding the critical energy density of matter corresponding to the big bounce is an 
open problem. 

It is tempting to apply fl50|) to the Planck scale. To make use of Planck's length Ipi :— 
y^HG/c 3 and Planck's energy density pp\ := c 5 /hG 2 , we multiply (|50|) by c 2 and recall that 
k 2 = 4vrG/3. Thus reads 

3 c 2 1 

Pmax ~ 8ttG 7 2 A 2 ' [ ' 

Substituting A = Ipi into (1BT1) gives p max = 3/87T7 2 p Pi ~ 2,07 pp\. Resolving (I5"T|) in terms 
of A makes possible finding A corresponding to pp\. We get A = a/3/87T7 2 Zp/ ~ 1,44 lp\. 



(We have used 7 ~ 0.24 determined in black hole entropy calculations |26l.l27||.) Surprisingly, 
the classical expression (15T!) fits the Planck scale. 

A natural next step is an examination of the energy scale of the big bounce at the quantum 



level. Preliminary calculations, [ri|, suggest that the energy scale would be described by 
the classical expression (I5TI) . 

B. Geometrical operators 

In the case the volume V = a 3 Vo of an elementary cell V is a cube, geometrical operators 
have simple dependence on canonical variables. Since the volume operator, V, is given by 

V(<j>) = v, (52) 

the area operator, A, reads 

A{4>) = v 2 '\ (53) 
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and the length operator, L, is found to be 



L{4>) = v l '\ 



(54) 



As far as we know, the geometrical operators have been considered so far only in the 



kinematical Hilbert space of the loop quantum gravity [28|, |29|, |30|, |3l|, |32| . We propose the 



examination of these operators in the physical Hilbert space of the loop quantum cosmology. 
It results from (1321) that we have 



K7A 



(55) 



so the geometrical operators are bounded from below by zero (as \\p<f,\ — > 0). 

An examination of the spectra of the geometrical operators at the quantum level is the 
next natural step. An interesting question is: Do these operators have the nonzero minimum 
and discrete eigenvalues? It is expected that the story will turn out to be similar to the 
case of the linear harmonic oscillator, where we have the nonzero ground-state energy and 
discrete energy levels. We present an answer to this intriguing question in our forthcoming 
paper lfl . 



VI. CONCLUSIONS 

We have shown that the resolution of the initial singularity of the flat FRW model with 
massless scalar field is due to the modification of the model at the classical level by making 
use of the loop geometry. The modification is parametrized by the continuous parameter 
A. Each value of A specifies the critical energy density of the scalar field corresponding to 
the big bounce. As there is no specific choice of A, the BB may occur at any low and high 
densities. The former case (big A) contradicts the data of observational cosmology (there 
was no BB in the near past!) and leads to weakly controlled modification of the expression 
for the curvature F£ b , i.e. gravitational part of the Hamiltonian (see the appendix). The 
latter case (small A) gives much better approximation for the classical Hamiltonian (see the 
appendix), but may easily lead to densities much higher than the Planck scale density, where 
the classical formalism is believed to be inadequate. Finding specific value of the parameter 
A, i.e. the energy scale specific to BB is an open problem. 

Our approach is quite different from the so-called effective or polymerization method (see, 
e.g. [33]), where the replacement (3 —>■ sin(A/3)/A in the Hamiltonian finishes the procedure 
of quantization. In our method this replacement has been done entirely at the classical 
level: Eq. (10) results from using an explicit form of the holonomy (5) in (8). Quantization 
consists in finding a self-adjoint representation of observables on the physical phase space 
and an examination of the spectra of these observables (ill . 12l . 



The elementary observables 0\ and O2 constitute a complete set of constants of motion 
in the constraint surface. They are used to parametrize the physical phase space and are 
"building blocks" for the compound observables like the energy density of the scalar field 
and the geometrical operators. So they have deep physical meaning. The role of 0\ and 
O2 becomes even more important at the quantum level as they enable finding quantum 



operators corresponding to the classical compound observables [13 . 

Our theoretical framework may be used for examination of the discreteness aspects of 
geometrical quantum operators, which may help in the selection of A. An extension of our 
formalism to the quantum level is straightforward. The algebra of observables is defined 
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in the physical phase space (hyper-surface in the kinematical phase space determined by 
the constraint equation). The carrier space of the self-adjoint representation of the algebra 
defines the physical Hilbert space. Examination of the eigenvalue problem of the length, 
area and volume operators (for fixed value of an evolution parameter they are observables) 
may lead to the specification of an unique A (or an interval of allowed values). Our next 
paper 16[ is devoted to examination of these problems. 

It may happen, however, that the value of the parameter A cannot be determined, for 
some reason, theoretically. The story may turn out to be similar to the case of the short- 
range repulsive part of the potential of the nucleon-nucleon interaction introduced to explain 
the scattering data [34] and the nuclear matter saturation of energy 35]. In such a case A will 



become a phenomenological variable parameterizing our ignorance of microscopic properties 
of the universe. 

An independent source of information on discreteness aspects of geometry is the obser- 
vational cosmology. The cosmic projects for the detection of gamma ray bursts may reveal 
that the velocity of cosmic photons depend on their wave lengths, which may be ascribed 
to the foamy nature of spacetime j36|, H3, Hl|. The detection of the primordial gravita- 
tional waves created at the big bounce may bring valuable information on the geometry of 
this phase (391 . 40, 41, 42|. The observational cosmology data may help to determine the 
phenomenological value of the parameter A. 
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APPENDIX A: HOLONOMY CORRECTIONS 

The curvature of SU{2) connection F k b = d a A b — dbA k + e^A l a A 3 b1 entering the expression 
([2]) for the gravitational part of the Hamiltonian, can be expressed in terms of holonomies. 
Using the mean-value and Stokes' theorems we have 

n F* h {x) « -i- / r k F* dx c A dx d ~ -^r (V exp ( f r k A k c alx c ) - l) , (Al) 

S ab Jct S ab V Jdu J 

where do is the boundary of a small surface a with center at x, and where s^ b := J dx a Adx b . 
The expression for F k b is exact but in the limit when we shrink the area enclosed by the loop 
da to zero. If we choose da in the form of the square with sides length /x, the expression 
for a small value of /i = /i has the form [l8| 

F*M = lim { - 2 Tr (^^J r k X °<4 + (A2) 

and we have 

F k ab = lim F k M. (A3) 

ii -^0 
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In the standard LQC the C(/U 4 ) holonomy corrections are ignored (see, e.g. 0, HI). 
It was found in 18|, ]l9|] that including higher order corrections leads to new curvature 
singularities different from the initial singularity and increases an ambiguity problem of 
loop cosmology. However, the holonomy corrections do not change the result that the big 
bounce is a consequence of the loopy nature of geometry (20| . 

Taking only the first term of (1A2j) leads to the simplest modification of gravity, but may 
be insufficient for the description of the inflationary phase. The choice of //o based on the 
expectation that Big Bounce should occur at the Planck scale Q has little justification (Io| . 
The significance of Planck's scale for quantum gravity seems to be rather a belief than proved 
result (see, e.g. j2l|). Heuristic reasoning playing game at the same time with Heisenberg's 
uncertainty principle, Schwarzschild's radius and process of measurement cannot replace a 
proof (see, e.g. [221]) . 



APPENDIX B: SYMPLECTIC FORM 

The symplectic form on the physical phase space Q may be obtained from the symplectic 
form on the kinematical phase space u by taking into account the constraint (1321 . 
The symplectic form corresponding to ffl3|) reads 



u = d<p A dps H \— d/3 Adv = d(j)A dp 6 H \— d/3 A dv. (Bl) 

47rG7 3kt7 



Making use of (132]) we get 



dv = KlX dp A - A ctg(A/3) d(3. (B2) 



Insertion of ( 1B2I) into ( 1B1I) gives 



Q = d<P A dp 6 + SSn } P ^ . d(3 A dps. (B3) 

3k sm(Ap) 



Since 

A <iarth(cos 
sin(A/3) ~ dp ' 

we have 



(B4) 



„ /,, sgn(ps) d&rth(cos\3) ,A , 

^ = # - Q ^ — — d/3 A rf^. B5 

\ 3k dp J 

On the other hand, due to (140]) . we have 

0i =p,, O 2 = 0-^M arth(cos(A/3)) . (B6 ) 

3k 

Therefore, 

= d0 2 A rfd. (B7) 
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Thus, the physical phase space may be parametrized by the variables 0\ and O2, and the 
corresponding Poisson bracket is given by (jUJ). 
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